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We study axion cosmology in a 5D Universe, in the case of fiat and warped extra dimension. The 
comparison between theoretical predictions and observations constrains the 5D axion decay constant 
and the 5D Planck mass, which has to be taken into account in building 5D axion models. The 
framework developed in this paper can be readily applied to other bulk fields in brane universes. 
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O ' I- INTRODUCTION 

(N ' 

One of the open issues of the Standard Model of particle physics is the so called "strong CP problem" [2, H, H, 0] . 
Because of the non-triviality of the QCD vacuum, the effective lagrangian must include the CP violating term 

»^G 2 . (1) 



Here g is the QCD gauge coupling constant, G 2 = G" G^ , where G" and G" are respectively the gluonic field 
strength and its dual, and 9 is a dimensionless free parameter. Measurements of the neutron electric dipole moment 
demand < 10~ 9 , which signals a very weak CP violation in the QCD sector, and immediately poses the question 
O i' of why can be so unnaturally small. At present, the most appealing solution seems to be the Peccei-Quinn (PQ) 
mechanism 0, @, 0, H, where the PQ symmetry, U[X)pq, is introduced and then spontaneously broken at low 
temperatures. That generates a Goldstone boson, the axion, which is capable of dynamically cancelling the effective 
\ 0-term and thus provides an explanation for the smallness of the 9 parameter. 
In the standard 4D theory, the axion anomalous coupling to gluons is 
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and the relation between the zero temperature axion mass, m a , and the axion decay constant, / a , is 
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where = 92 MeV and — 135 MeV are respectively the pion decay constant and mass, while z = m u /md ss 0.6 
J^j ' is the up to down quark mass ratio. 

The bounds on the axion decay constant depend on the exact axion model, but the allowed range is essentially [T(| 

10 8 GeV < f a < 10 12 GeV . (4) 

The lower bound can be deduced from considerations on stellar cooling or on the neutrino burst of the supernova 
SN1987A, while the upper limit from the requirement that the present axion relic density does not exceed the dark 
matter energy density. To be more precise, supernovae forbid only the range 10 5 GeV < f a < 10 8 GeV, but from other 
considerations (e.g. the overproduction of axions in the early Universe by thermal processes) we find f a > 10 5 GeV. 
Using eq. ([3]), the constraint ([4]) can be written as 

1CT 5 eV < m a < 1CT 1 eV . (5) 

Axions can be also produced non-thermally via the following (misalignment) mechanism [IT| . At high temperatures, 
the cosmic energy density of the axion field is p a = m a (T) a 2 . Here m a (T) is the thermal axion mass, which has to 
be calculated by considering QCD instanton effects at high temperature. For T w 1 GeV, one finds 



»ori = n.lm„(^j , (6) 
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where a 2 ~ (0/ a ) 2 is the value of the axion field. As the temperature of the Universe decreases, the axion mass 
increases and when it becomes of order the Hubble parameter, H = 1.7 y/g*T /Mpi, the axion field starts oscillating 
and converts its energy into axions. This happens at the temperature T m , given by m a (T m ) w 3H(T m ). The axion 
number density to entropy density ratio at T m is 

= ra a (T m )9f„/2 
s s 

where s = (27r 2 /45) g*T^ (<?» counts the relativistic degrees of freedom at a given T). The present value of £l a = 
Pa/ Pc can be obtained by multiplying eq. ((7|) by the zero temperature axion mass and the present entropy density 
so = 2970 cm -3 and by dividing over the critical energy density p c = 2 • 10 -29 h 2 g cm -3 . The result is 

n a fc 2 «o.6e? n (y^v) > ( g ) 

and implies f a < 10 12 GeV for Qi n ~ 1 (Qi n is expressed in radians). Actually @i n may also be much smaller, 
say at the level of 10~ 3 , in order to have f a close to the GUT scale. However, if the PQ phase transition occurred 
before inflation, @i n is limited by the bound from isocurvature perturbation on the CMB 112 . HH, [l4| • But the bound 
depends on the field value of PQ scalar during inflation which may be different from f a Thus the isocurvature 
constraint is model-dependent and we do not consider this constraint in the present paper. 

The scope of this article is to study the cosmology of a 5D axion field in a 5D universe, in which our Universe is 
thought of an embedded 4D hypersurface. In particular, we will focus our attention on 5D braneworld models, which 
are the most studied and discussed in the literature and which represent the simplest extension to our 4D Universe. 
In the next section, we review the basic ingredients necessary for our analysis, that is, basic facts of 5D axion models 
and the master equations governing the evolution of the 5D universe. In section IIII1 we study the most important 
axion production mechanisms (thermal inelastic scattering, time varying metric, coherent oscillation of the axion field) 
and show for what values of the free parameters of the theory each process turns out to be the most efficient. In 
section IPVl we present the constraints we can deduce from astrophysical and cosmological arguments in our picture. 
Summary and conclusions are reported in section [Vl 

Conventions: Throughout the paper we use natural units h = c = ks = 1. As for the notation for the different 
gravity scales, M Pi = 1.2 • 10 16 TeV is the 4D Planck mass, M 4 = 2.4 • 10 15 TeV the reduced 4D Planck mass and Af 5 
the reduced 5D Planck mass. 



II. 5D AXION MODELS 

In a 5D spacetime, where the Standard Model (SM) particles are confined to the 4D brane (see [H|), while gravitons 
and axions live in the bulk, one expects a much richer phenomenology compared to the 4D world [171 [18[ , mainly due 
to: 

• Now the "fundamental PQ scale" is the 5D symmetry breaking scale f$. The effective 4D scale, /4, depends on 
the former and on the shape and the size of the extra dimension; 

• From our low energy 4D viewpoint, a bulk field is seen as discrete or continuous tower of Kaluza-Klein (KK) 
states; 

• The mass of the zero mode axion may not depend on the PQ scale in the usual way m a ~ f-Kin^/ f a . To be 
more precise, in the toy-models of refs. [13, Ell], one finds that the mass of the mode is roughly min(l/i?, m a ), 
where R is the size of the extra dimension; 

• If the highest temperature at which the Universe attained thermal equilibrium is higher than the so-called T t 
(see below), the Universe experiences a period of braneworld expansion 0, [2(| l2lj . This, for example, could 
imply an efficient gravitational production of axions. Notice that, in fact, there need not be equilibrium thermal 
plasma for a non-standard expansion epoch to occur; however, for the phenomenology we are studying here, 
this is the most relevant case. 

Some of these considerations apply, of course, to any bulk field, and can be easily extended to other 5D models. 
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Flat extra dimension 



A 5D axion in a flat compact fifth dimension (a circle) is discussed in ref. [l7| (note that here we adopt a slightly 
different notation). The 5D axion action is 



2 y ' 32tt 2 ^3/2 



(9) 



Saxion — d x dy 

where a is the 5D axion field, whose energy dimension is 3/2. Expanding a into Fourier modes 

a(x, y) = J— V a n {x) b n {y) (10) 
V2ttR n 

and integrating over the extra dimension, we get the effective 4D lagrangian. The effective 4D PQ scale is = 
(27ri?/5) 1 / 2 /5 and, for Rf§ 1, it can be much higher than the 5D decay constant f$ (the picture is like that of the 
4D and 5D Planck masses, which are related to each other by M4 = {2ttRM^) 1 ^ 2 M5). 

Physical phenomena can be easily described and interpreted using a 4D or a 5D point of view. In the first case, 
we have KK axions with coupling constant f± and mass splitting 1/R (this is valid for high KK modes, and is 
only an approximation for lower states, see below): the axion production rate at the energy E is given by the the 
production rate for a single state, ~ E 3 /f%, times the number of states with mass smaller than E, i.e. RE if RE 3> 1, 

r ~ E^ • Adopting the 5D point of view, there is just one axion field, but the axion coupling is I//5 2 - The rate 
is T ~ E //f, and the two descriptions are equivalent, since /| ~ Rfi- 

Let's write down the few equations that will enable us to switch from one set of parameters to another (in TcV 
units). 

1/R ~ 1.1 • 10 -30 Mf ; (11) 

h - 2.4-10 15 (/ 5 /M 5 ) 3 / 2 ; (12) 

m a (eV) ~ 2.5- 10- 12 (M 5 // 5 ) 3/2 ; (13) 

mo = min{l/i?, m a } ■ (14) 

Finally, in general the mass of each KK mode is given approximately by m n = mg + n/R, where R^ 1 = 2ttM^/M^ 
is the inverse of the size of one compact extra dimension. This expression of course is not universally true, since an 
effective potential for the 5D field may appear upon dimensional reduction, and therefore contribute to the effective 
mass of the KK mode. This is indeed what happens in the extra dimensional axion models at hand, but nevertheless, 
it is not difficult to convince oneself that at high n this parametrisation provides a very good approximation. This 
applies to the warped case as well. 



B. Warped extra dimension 

The case of one warped extra dimension is discussed in ref. [18j |. However, the picture is quite similar to the flat 
case and likely there are no relevant differences. 

If the extra dimension is warped, the structure of the KK tower turns out to be subtly different. In the case of two 
branes, typically one finds that the mass of the KK mode n is 

m n = m + kxneT^^ , (15) 

where k = (1 — exp(— 2irkR))M§ /M% is the AdS^ curvature, x n is the n-th root of the first order Bessel function J\ 
and ttR is the size of the orbifold. The mass splitting is Am ~ 3fc exp(— -KRk), and for later convenience we define 
the function F(kR) = (1 - exp(-27rfci?))/ exp(7rfci?). 

Thus, some care should be taken in working out the formulas for the abundances, as the expressions for the mass 
gaps are different. However, as it will be shown below, in the end the picture is almost unchanged. 



C. Brane cosmology 

As it has been noticed at the beginning of this section, one of the main differences between 4D and 5D cosmologies 
lies on the different relation linking time and temperature. Indeed, if the SM lived on a four dimensional Friedman- 
Robertson- Walker hypersurface (the brane), embedded in an extra dimensional spacetime, the early universe would 



4 



admit an epoch of non-standard expansion [2l|. Several such models have been built in the last few years, see [22|, 
23, 24] for sonic of the original proposals. These models show a peculiar feature when their cosmology is investigated. 
Indeed, a general feature of the 5D embedding is that on each 4D brane the Friedman equation has the form (l9l. [20. l2ll] : 

3M| V r 2A 



* = &[l + &), (16) 



where A is the tension of the brane, which is related to the five dimensional Planck mass as A = 6M5/M4. This 
equation says that at high energy densities the expansion of the universe was much faster than at later times, and 
went as T 4 instead of T 2 , together with the unknown parameter M5: the smaller M5, the faster the expansion. Here 
we have ignored further terms in the effective 4D equation, but they are required to be very small by observation, see 
for instance [25| . 

At this point it is convenient to define a "transition" temperature T t from standard cosmology to brane one, which 
can be extracted from p — 2A: 

360 \ 1/2 M? 



T t ?= -5— TTi (17) 



Ma ' 



where g* = g*(T) counts the relativistic degrees of freedom at a given temperature T. If the dominant component 
of the universe is not radiation, then this "temperature" approximately means the fourth root of the energy density, 
and parametrises the epoch at which the transition occurs. 

One basic requirement is that the Universe expands as H 2 oc p at the BBN, which means T t > 1 MeV, or 
M 5 > 10 TeV. 

III. PRODUCTION MECHANISMS 
A. Thermal Production 

In this section we will be analysing the conditions under which a thermal population of (KK) axions is born in the 
early Universe. The results obtained here are not only relevant for the thermally generated abundances, but also in 
dealing with other non-thermal production mechanisms (gravitational particle production, misplacement mechanism, 
etc), as non-thermal abundances generated before the epoch at which axions thermalise will be phagocytosed by the 
thermal bath. 

The most important process which thermalises axions in the early Universe (assuming inflation takes place at higher 
scales) is their interaction with QCD matter. The rate of these processes (see (26|) can be estimated as 

r ~ 7 • 10~ 6 T 3 //J ~ 1.2 ■ 10~ 36 T 3 A/|// 5 3 , (18) 

and the expansion rates of the Universe are given by (in TeV units) 

H std ~ 1.4 • 1{)' W ^T 2 , (19) 
H bc ~ 6 • 10- 2 .g,T7M 3 . (20) 

It is straightforward then to realise that thermal equilibrium took place when T > H, which implies that 

T atd % T s°t d - 1-1 ' l0 W V9ji/Ml , (21) 
T e b q c <T b D c ^2- IQ-^g^Ml/fi , (22) 

where T® d and are the decoupling temperatures, respectively in the standard and braneworld expansion period. 
The transition temperature between the two different regimes (from H 2 ~ p 2 at high energies to H 2 ~ p at late times) 
happens at 

T t ~ 5 • lQ- 8 gZ 1/4 M^ /2 . (23) 
In order for the two equations (|2"Tj) and (f2"2"|) to make sense, the following condition must be satisfied: 

T b D c >T t >T° d => f 5 <S-10- m g- 1/4 M 3 5 /2 , (24) 
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or the band of temperature for which thermal equilibrium is realised would shrink to nothing. Notice further that if 
the reheating temperature or the 5D PQ breaking scale were smaller than the upper limit T^?, then they would be 
the actual upper (in temperature) limit at which thermal equilibrium ceases to be realised. Of course, in the case the 
reheating temperature is even smaller than T® d , axions were never in equilibrium after inflation. 
Now, for each mode in thermal equilibrium which decouples when still relativistic, we have: 

Tl 

Y n = — = 0.278/ 5 » . (25) 
s 

Here, as previously, <?, is the number of relativistic degrees of freedom in equilibrium with the thermal plasma. In 
our case (formula (|25[> ) this is taken at the decoupling temperature. Moreover, if there are many KK axion states 
in thermal equilibrium, this number is given by gsu + tikk with tikk accounting only for the KK states which are 
relativistic and in equilibrium, and gsM = 106.75 for the SM, or about twice as much for its supersymmetric version. 

In order to render the analysis traceable and clear, it is most useful at this point to specify a working value for f^. 
With a little foresight, we pick the value / 4 = 10 7 TeV: we will comment on the dependence of our results on this 
choice every time it will turn out to be relevant. 

By fixing a value for f± we automatically tie /g to M5 as 

h - 3 • 10- 6 M 5 . (26) 

Flat fifth dimension - Since from this point on we will need the explicit expression for the mass gap, it is sensible 
to split the analysis in two, flat and warped 5D axion models. 

First of all, let us look now at the behaviour of the effective number of light degrees of freedom with the 5D Planck 
mass. As we are ultimately interested in this number at decoupling, we write the number of relativistic axion modes at 
decoupling as n^ K — T^ d R. There is therefore a critical value for M5 at which the equilibrium KK axions overwhelm 
the SM degrees of freedom; this happens around M5 w 5.6 • 10 10 TeV. This special value of M5 marks the point at 
which, at decoupling, the number of equilibrium KK axions reaches (or drops below) the number of SM degrees of 
freedom. Consequently, the decoupling temperature T^ d becomes 

Tf td ~ 3 • 10 4 TeV when M 5 > 6 • 10 10 TeV , (27) 
Tf td ~ 4 • 10 36 A/ 5 ~ 3 TeV when 1.2 • 10 10 TeV < M 5 < 6 • 10 10 TeV , (28) 

where the first line applies to ~ gsM and the second case is instead g„ ~ T® d R. Hence, the number of equilibrium 
KK states that decouple at T° d is almost always smaller than the SM ones, unless 1.2 • 10 10 TeV < M 5 < 6 • 10 10 TeV. 

The picture thus is as follows. In the early Universe, even before T t we have a population of thermal axions because, 
as long as (|24p holds, the window for thermalised KK axions is open, and goes from T® c down to T^ d , crossing T t , The 
number of states that are still relativistic and in thermal equilibrium at the lowest decoupling temperature depends on 
the value of M5 as shown in (|27[) and (|28|) . Since we have assumed that the cross section is independent on the mass 
of the KK state, then all these states will come out of equilibrium at the same time when T = T® d , and all of them 
will be by definition relativistic then. Of course also axions of higher mass were in thermal equilibrium before the 
collective decoupling, but they were non relativistic and their abundances are therefore irrelevant due to Boltzmann 
exponential suppression. Once the axions decouple, their yield variable remains pretty much the same until today, 
unless some modes decay (or other production mechanisms) turn on. 

Given these results and this picture in mind, we can rewrite the yield variable at decoupling as 

{3 • 10~ 3 when M 5 > 6 • 10 10 TeV 

(29) 
9 • 10~ 68 Ml when 1.2 • 10 10 TeV < M 5 < 6 • 10 10 TeV 

which is valid per each d.o.f. up to tlkk- 

Probably the easiest way to analyse the impact of KK towers in late time cosmology is to split the tower into mass 
bands according to the constraints one is going to look at. This is so because different constraints corresponds to 
different decay time ranges, which in turn, once the zeroth mode decay width has been given, depend almost solely 
on the masses. Once the interesting lower and upper limits have been identified, one is in a position to scrutiny the 
number of modes in that given range, and investigate their total impact on a particular observable. Moreover, it is 
reasonable to expect that only large number of states for each band could provide significantly different constraints 
compared to an ordinary 4D particle with mass within the given band. This last consideration, combined with the 
fact that the zero mode is supposed to be light, implies that the highest KK mode N in the band will satisfy JV>^o, 
with Nq the lightest mode in the same band. This scheme was endorsed in [27], HH to investigate thermal production 
of 5D gravitinos, and in [2^| for gravitationally produced particles (see also below). 

By looking at (|29p once can readily evince that single modes, in order to be safe, demand: 
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• Small mass, therefore a small M5 

• A high T® d which implies more dilution afterwards. This wants a high /g (not displayed in our formulas), 
although there is a threshold beyond which thermal equilibrium is never realised; it also requires a small M5, 
again. 

• Much dilution, that is, the largest Ukk available, which again supports small M5 and high f$. 

The situation changes drastically when a summation is involved. Indeed in this case the solution to overproduction 
of axions is a widening of the mass gap between different states, which implies on the one hand less dilution (which is 
bad) , but on the other hand precludes an enormous number of states to be summed over and contribute to the total 
energy density today. This can be put in formulas as 

M 5 > 6 • 10 10 TeV 

(30) 

1.2 • 10 10 TeV < M 5 < 6 • 10 10 TeV 

which (remember that ./V oc M^~ 3 ) shows the inverse dependence on M5, thereby imposing a lower limit on M5. Notice 
that this is not the case when the number of equilibrium states is dominated by KK axions, but it will be true for the 
other production mechanisms. 

Warped fifth dimension - If the extra dimension is warped the mass gap is much different than in a flat bulk, 
difference which is traceable by following the effects of the function F(kR) in the previous formulas. 

In order to have an idea of what happens in this case, let's work out the equilibrium conditions (analogously to 
eq. (J2U)). Firstly, it is necessary to know which degrees of freedom dominate at decoupling. In doing so, we separate 
the problem in two, by noticing that while the highest value F can take (that is, around 0.4) would leave the coefficients 
unchanged, for much smaller F's (equivalent to more pronounced warping) the situation is dramatically different. 

Indeed, while in the first case there are only tiny numerical differences, if the warping factor grows there are more 
and more KK axions in thermal equilibrium (for fixed M5) and if we take kR — 11 the KK equilibrium d.o.f. always 
dominate the SM ones. In this second case, thermal equilibrium is realised as long as M5 > 8.1 • 10 13 TeV, and the 
predicted thermal abundance is 

Y n {t D )~ 1.9 -10- 98 M 5 6 . (31) 

This, as long as M5 > 10 13 TeV, would be a disastrous amount of KK particles (> 10~ 14 per each mode), but they 
luckily decay just before the onset of the BBN and, therefore, are reasonably harmless. 

B. Gravitational production 

Gravitational particle production in an ordinarily expanding Friedman-Robertson- Walker (FRW) Universe domi- 
nated by dust-like or radiative matter is known to be a very poorly efficient particle creation mechanism [30L l3l| . 
However, as it has been recently pointed out in ref. (32j |. if the four dimensional FRW Universe were to be embedded 
into a higher dimensional spacetime, then even dust-like or radiation dominated FRW Universes would be able to 
inject gravitationally produced particles into the plasma in a sizeable way. This is easily understood as a consequence 
of the possibility that the actual scale of - 5D - gravitational interaction be tuned to much lower values, thereby 
drastically enhancing the coupling with matter. 

In what follows we simply borrow the results of [29l |32| , referring the interested reader to these papers for the 
details. 

Flat fifth dimension - Having obtained explicit formulas for the KK masses, one can readily write down the 
gravitationally produced abundance associated to each KK mode. In the limit for which the zero mode is light 
compared to n/R, the yield variable turns out to be 

Y n ~ 1.2- KT 71 ™ 9 / 4 M 5 9/2 . (32) 

In order for all this to be consistent, we need to ensure that 

T f >T t =*> n > 0.3 , (33) 

What this means (recall it is an order of magnitude estimate) is that the result ([3^]) does not hold for the first handful 
of modes, as the quantities fl n and Y n have been deduced assuming H = p/6Af 3 , while for small n we would find 



x I 1.4 • 10 -33 iV 2 Mf TeV when 



•S 



n=0 



5 • 10" 98 V 2 Ml TeV when 
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Tf < T t . For these light modes, the gravitational production stops at the transition temperature T t or, if very light, 
during the standard expansion, thereby rendering their abundances negligible. 

With the single mode contributions at hand, it is now straightforward to obtain the overall yield variable for each 
band by just summing over the modes (up to the highest mode N) and then discarding the (smaller) contribution of 
the lowest state, the result of which being 

^ = Y m n Y n ~ 3 • 1(T 102 7V 17 / 4 M} 5/2 . (34) 

n 

Warped fifth dimension - The main difference arising in this case is the different expression for the KK masses, 
as the Friedman equation is basically the same of the case with a flat extra dimension and we can still use the modified 
Friedman equation. It is thus fairly straightforward to repeat the steps undertaken for the flat case, with result 

y„~2-10- 72 n 9 / 4 F 9 / 4 Af 5 9/2 , (35) 

in place of eq. (13"2"1) . 

Similarly to the case of thermal axions, the scenario is more or less unchanged in this picture, unless one employs 
very small values for F(kR). Indeed this function is always smaller than ~ 0.4, for which value the abundances and 
consistency constraints turn out to be very similar to those previously obtained for a flat extra dimension. If one 
instead chooses to work with much smaller F(kR), such as for kR = 11, then the mass splitting becomes extremely 
tiny, unless the five dimensional Planck mass is pushed all the way up to the 4D M4. Therefore, the mass gap drops 
by a factor of 2/F, and the abundances of produced particles increase by a huge factor 1/F. The compensation in M 5 
is proportional to F~ 4 ^ 21 which, again for kR = 11, is around 700, although in highly warped models it is customary 
to safely (as far as this mechanism is concerned) take M5 ~ M4. 

Along these lines one can now transfer to the warped solution the results which will be presented below for the flat 
scenario. 



C. The misalignment mechanism 

Flat fifth dimension - In the standard scenario, (often) the most efficient production mechanism of axions in the 
early Universe is associated to field oscillations. Such a mechanism works because massless (or light, that is H S> m) 
scalars are unstable in dS spacetime and during inflation can acquire a large vacuum expectation value (vev). Then, 
when the expansion rate of the Universe becomes smaller than m, the field goes to the minimum of its potential and 
converts its energy into particles. In order to know how many modes have acquired a large vev, we should know the 
value of the Hubble parameter during inflation, Hj. The number of modes would be N w RHj. 

Other important inputs are the temperature dependent axion masses (temperature corrections are likely significant 
for the zero and the lightest KK modes, while are more probably irrelevant for the other ones) and the kind of 
expansion (standard or brane regime?). In other words, the estimate of the axion abundance depends on several 
parameters, which we will try to track in our analysis. 

Let us assume that the PQ symmetry is already broken after inflation. In this case, the axion field is homogeneous 
over large distances and we have to consider only the zero momentum mode 1 . 

The abundance of KK modes can be evaluated in a very similar way to that described for the zero mode, but now 
using the KK mass m n = n/R = 2imM§/M%, and assuming that the thermal mass is not relevant (which will be the 
case for heavy KK states, which in turn are the states we are mostly dealing with). Since we are interested only in 
the heavy KK modes, we will not go into the details of the computation of the thermal mass. For > T t , the n-th 
KK mode oscillates during the braneworld period and we get 

I£ - 1.6- 10" 8 n 1/4 M 3/2 TeV, (36) 
Y n = 6 • 10V/ 4 e 2 „ A/ 5 " 3/2 . (37) 
For < T t , the n-th mode oscillates when the Universe expands in the standard way and we find 

T™ = 1.6- 10" 8 n 1/2 M 3/2 TeV, (38) 
Y n = 6- 10 6 n- 1/2 e 2 „M«r 3/2 . (39) 



If the PQ phase transition is after inflation, that is not true and we have to consider non-zero momentum modes as well, whose 
contribution to the axion cosmological energy density would be similar to the one of the zero momentum mode. So, neglecting such a 
possibility we find at most a more conservative bound. 
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The total axion abundance can be easily deduced by summing over all the KK modes which acquired a large vev during 
the inflation period. That is determined by the exact model of braneworld inflation. A more accurate estimate of 
axion production from coherent oscillation would require a numerical integration of the field equation in the expanding 
background. Nevertheless, such an effort would be likely useless, because the most important parameter is the initial 
misalignment angle Oi„ which is unknown. 

Notice also that the temperature T m is always (for our choice of f&) above the transition temperature, except for 
the first few KK modes. Therefore, in the discussion of the constraints, we will employ only the abundance calculated 
from ([571). 

A comment is in order here, concerning the analysis of relic axion oscillations presented in [ItJ ■ There a mechanism 
for enhancing the rate of dissipation of the total (zero mode plus KK states) energy density stored in the axion 
condensate is exploited for rendering the 5D axion model possibly even safer than the standard 4D scenario. 

This mechanisms relies on one fundamental assumption, that is that all KK modes sits at the minimum of the 
potential (<di n — 0) when the zero mode begins its oscillations. In this case the zero mode triggers these higher KK 
modes into oscillations, and therefore transfers some of its energy density into them. As is well known, heavier axions 
during the subsequent evolution of the Universe lose energy more efficiently than light ones, which is why axions 
produced by the misalignment mechanism put a bound on the invisibility of the axion (that is, require a minimum 
mass). Hence, the net result is that KK modes help the dissipation of the total energy density in this model. 

However, since heavy KK modes are expected to decay the dissipation happens through the first few KK states. 
But these states are also the states that are likely to have some initial displacement from the bottom of their potentials 
(which appears even before the QCD transition due to the mass matrix between KK states). This initial displacement 
is not accounted for in [l7j . and we therefore expect that the benefits of this mechanism will not hold in a more 
realistic case. 

One last brief note before moving to the warped scenario. In the analysis of [17| the effects of the brane regime 
were not included; indeed if one looks at the values of fi and Mg needed to have an enhanced dissipation rate, one 
finds that, for / 4 « 10 7 TeV, the 5D Planck mass (assuming only on extra dimension) needs to be below 500 TeV, 
but we know that then the zero mode will start oscillating before T t , and the equation of motion will therefore be 
different (due to the different dependence of H upon time). 

We can now write down the total abundance, keeping Oj„ as an unspecified parameter, as 

^ ~ 3 • 10- 24 N 9 ^e 2 ln M 3 5 /2 . (40) 

Warped fifth dimension - Like for the other production mechanisms, here we have just to replace the KK 
mass of the flat model, m n — 2nnM^ /M 2 , with the one of the warped model, m n — 3F(kR)nMj- jM\. The same 
considerations presented for the flat case apply to this scenario as well, and will change the results of ref. . 



IV. CONSTRAINTS 



In this section we will be analysing the bounds coming from different cosmological arguments, such as the density 
of Dark Matter particles today, BBN constraints, CMB distortion, X-rays, effective number of neutrinos, and recom- 
bination. In order to do that, we need to know the lifetime of the KK axions. For masses larger than about 100 MeV, 
the main decay mode is usually into hadrons and the lifetime is 

For smaller masses, the decay into hadrons is kinematically forbidden, and the main channel is a — > 27. The lifetime 
is given by eq. (|4Tj) with a s replaced by a em : 

— » 3 '0 3 (lo4v) % - (42) 

In the end of this section we will briefly reconsider the supernova bound as well. All results refer to the flat extra 
dimension case, but it is straightforward to export the results to the corresponding warped scenario, thanks to the 
explicit formulas for the abundances computed in the previous section. 

Before embarking on to the discussion of all these limits, a general remark on thermal abundances is in order. We 
have seen that thermal equilibrium for the axion (at fi = 10 7 TeV) occurs only if M5 > 10 10 TeV or so. However, this 
implies automatically that the first KK mode will have a mass of order TeV, which implies that it will decay much 
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faster than Is, and therefore provide no constraint at all. This will not in general be true for different values of 
and more in general for other bulk fields, but in this specific case thermal abundances of KK states are automatically 
made safe. 

A description of each way we can constrain the abundance of KK axions follows, together with the rough estimates 
for the limits on the 5D Planck mass. A summary table with the more precise limits is presented at the end. 

A. Cold relics 

Axions whose lifetime exceeds that of the Universe to ~ 10 17 -^10 18 s, will contribute to the energy density of invisible 
matter today, which is constrained by observations to be not in excess of Vl a h 2 w 0.12. This can be translated as 

^<5.10- 13 TeV(^M^l . (43 ) 
s ~ V °- 12 / 

The number of KK states which have survived long enough to be counted in this quantity are all the light ones up 
to N ~ 3 • 10 22 Mr" 3 , and their contributions due to gravitational and misalignment abundances are 

& ~ 1.5 • 10- 6 M 5 - 21/4 , (44) 
^ ~ 1.3-10 27 e 2 „M- 21/4 , (45) 

which in turn means that the 5D Planck mass has to be bigger than about 17 TeV and 3 • 10 7 0f„ 1 TeV, respectively. 

Notice that these constraints, as all those that follow, are valid only if there effectively is a densely populated 
tower of states within the mass range considered, fact which can be easily checked by comparing the highest mass 
for each band with the minimum mass gap required by the corresponding constraint. As an example (which can be 
reiterate in the following discussion), in the case just worked out the highest mass involved is of order 40 keV, and 

the minimum mass gaps are around 10 -15 eV and 40 O^ 21 keV for gravitational abundances and misalignment ones, 
respectively. Therefore, while the result obtained from gravitational production is surely valid, the one referring to 
the misalignment mechanism is good only for small angles Qi n . This is a general result which is found to hold in the 
forthcoming sections as well. 

B. BBN 

If a particle lives long enough to witness the synthetisation of the light elements (Big Bang Nucleosynthesis, or 
BBN), but decays during this process or shortly thereafter, the entropy density injected into the equilibrium plasma 
can potentially alter the final abundances of these elements, thereby causing conflicts with the good matching between 
theoretical predictions and observations. Limits on the amount of energy density to entropy ratio which can be safely 
transferred to the plasma are usually among the strongest ones for particles whose lifetime is in the range of interest 
t 1 s to t 10 13 s. In working out the constraints on the KK axion abundances we closely follow the schematisation of 
ref. [33| . and we denote B 1 and £?had the branching ratios into photons and hadrons, respectively. 

The four possibilities we consider are: 1) = 1 with 10 4 s < r < 10 7 s; 2) B 7 = 1 with 10 7 s < r < 10 11 s; 3) 
Shad = 1 with Is < r < 10 4 s, and finally; 4) B had = 1 with 10 4 s < r < 10 n s. This last band will not be able to 
constrain the parameters of the extra dimensions, since the hadronic axions with masses above 100 MeV decay before 
10 4 s. Finally, notice that, although limits can be obtained for decays down to r 10 13 s, the existence of the 4.5 MeV 
threshold (below which there is not enough energy to dissociate D + 7 — > n + p) , effectively cuts the lower end of the 
bands around 10 11 s. 

Moreover, for the sake of simplicity, the allowed total p a /s have been flattened over the timescales (or mass ranges) 
pertaining to each band, which means that our constraints are actually conservative. 

The three possibilities are reported in table Q] where we can see that gravitational particle production limits the 
lower value of M5 around 10 4 TeV to 10 5 TeV, while, although some care should be taken when reading these limits, 
misalignment axions tend to push M5 up to values for which the first KK mode should have decayed before Is. 
However, if the misalignment angle is small enough (smaller than 10 -4 for example), then the strongest limits are 
comparable to those coming from gravitational production. 
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C. CMB distortions 



If the KK states have a lifetime between 10 6 s and 10 13 s, then their decay may affect the blackbody spectrum 
of CMB. Indeed, such decays would transfer energy to the plasma at the epoch when processes such as the double- 
Compton scattering, necessary to preserve the blackbody spectrum of the CMB, become inefficient. The distortion 
of the spectrum is characterised by the chemical potential n at times greater than 10 9 s when the energy transfer by 
the Compton scattering is efficient, and Compton y-parameter at later epoch. Constraints on these quantities can be 
read off [34j | . and, once translated into upper bounds on the energy density to entropy ratio, become 

— <3-10~ 14 TeV for 10 6 s < r < 10 9 s , (46) 
s 

^<7.10- 16 TeV for 10 9 s < r < 10 13 s . (47) 
s 

The resulting lower limits on the 5D gravity scale M5 are similar to those extracted from BBN: around 10 4 TeV 
for the gravitational abundances, and 10 9 ©^{ 21 TeV for the KK axions produced via the misalignment mechanism. 
Notice that, as seen in the summary table below, KK axions that decay around 10 6 s have masses slightly above the 
threshold for the decay into hadrons, so the limits we are quoting are valid only for i?had <C 1. 



D. Effective number of neutrinos 



The expansion rate of the Universe depends on its total energy density. The study of phenomena which are sensitive 
to the Hubble parameter can thus put constraints on the amount of new particles. The best example is the BBN, 
where the primordial production of light nuclei depends on the expansion rate when the Universe temperature was 
in the range 1 MeV — 10 keV, that is, at the time t » 1 — 10 4 s (see e.g. ref. [SI])- The bound that can be deduced 
is often expressed in the literature as the effective number of neutrinos and current data demand no more than one 
extra neutrino species: AiV„ < 1. 

The KK axions which are stable or quasi-stable at the BBN are the ones which can not decay into hadrons and 
whose mass is not larger than about 100 MeV. If we assume 1/R -C 100 MeV, i.e. M5 <C 5 • 10 s TeV, the highest 
KK mode is N ~ 10 26 M 5 -3 . The total KK axion energy density to entropy density ratio due to gravitational and 
misalignment productions are respectively 

^ ~ 7-10 8 M" 21/4 , (48) 
^ ~ 7-10 34 e2 ri Af- 21/4 . (49) 

Such a quantity must be smaller than \.3g* l {T pro d)TBBN , where g*(T pro d) counts the number of degrees of freedom 
at the time of axion production, while Tbbn is the temperature at the BBN. So, BBN demands M 5 > 1.52 ■ 10 3 TeV, 
if we consider the mechanism of gravitational production, and M5 > 1.4 • 10 8 TeV in the case of coherent oscillations 
of the axion field = 7r). 



E. X— ray background 

Photons with energy in the range 1 keV — 1 TeV produced at time t > 10 13 s can contribute to the diffuse X-ray 
background. The flux of photons from the two photon decay of a particle of mass m and lifetime r is found to be [36j 



*° , B^n(z) dE',, 
dt' ,7 y k -7T7 2(5 (E' - m/2) 



n 



B~ ( s Y 
2tttH, 



i) 



t(1 + z) 3 dE 



f(m/2E) exp 



3F rOi /2 



In 



f(m/2E)fl 1 / z - 1 tt 1 / 2 + 1 



1 f^ /2 ll 



(50) 



where B 1 is the branching ratio into 27, n(z) is the particle number density at redshift z, E' is the photon energy at 
the instant of production, E is the present photon energy, 



f(m/2E) = 



/ vn \ ~- 



-1/2 



(51) 
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and in the last step we assumed that the Universe is flat, i.e. Q\ + f2 m = 1. The estimated photon flux must be 
smaller than the observed one [37], [H[ : 

/ E \ ~°' 4 

$obs = 8 (lkev) cm^s^sr- 1 , (52) 



in the energy range E m 0.2 — 25 keV, and 



E 



$ofcs = 6 ' 10 ^ I ^McV ) Cm " 2 Sr_1 ' (53) 



for E w 25 keV - 4 MeV. 

The strongest constraint comes from KK axions with lifetime much shorter than the present age of the Universe. 
In this case, we can assume that all the axions of the level n decayed at the time t — r„ and we find that the photons 
produced in the decay have today energy 

For example, photons of the X-ray background with energy in the range 1 — 2 keV might have be produced by the 
decay of cosmological KK axions with mass 250 — 500 keV. Since the amount of axions produced via gravitational or 
misalignment mechanism are respectively 



•s 



0.1M~ 21/4 , (55) 



^ ~ 5-10 29 e 2 „Af~ 21/4 . (56) 



Requiring 

< 



Fobs (E) dE , (57) 

1 keV 



s ~ 3 so 

we get the constraint M5 > 7 • 10 3 TeV from gravitationally produced axions. 

For the misalignment mechanism, if we use equation (|5T|) , we find M5 > 5 • 1O 9 0*,( 21 TeV. This result however is 
valid only for small angles, as when Qi n approaches it the mass of the first KK state reaches (and surpasses) 500 keV, 
which means that today the photons produced in the decay are below 1 keV (see table H}. 

F. Reionisation 

Photons produced by particles decayed after recombination (r > 10 13 s) are redshifted due to the expansion of the 
Universe and may thus leave the transparency window £sl keV — 1 TeV. If this were the case, they would interact 
with the intergalactic medium and provide an extra source for reionisation. The abundance of the parent particles 
can be constrained by demanding that the contribution from these photons is still consistent with the optical depth 
of the last scattering surface obtained from WMAP. 

In our case, reionisation can constrain the abundance of KK states with mass around 1 MeV, which decayed at 
t 10 13 s into two photons, whose energy today is expected to be about 0.5 keV. Assuming 1/R <C 1 MeV, the 
fraction of energy density of KK axions with mass around 1 MeV is 

il g w 1.6- 10 13 A/~ 21/4 , 

n m a 1.2.10 42 e 2 n M- 21/4 , (58) 
respectively from gravitational and misalignment production. Observational data require that [39j 2 

B 7 n g , m < 10- 11 . (59) 

So, assuming £> 7 w 1, we find respectively M 5 > 3 • 10 4 TeV and M 5 > 1.1 • 10 10 6- ? { 21 TeV. As already pointed out, 
the last bound is only good for small initial misalignment, while large angles imply that the first KK state would need 
to be heavier than about 1 MeV. 



2 The parameter f of ref. [39ll can be read as fig^/n^m here. 
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Gravitational production 





Abundance (times Mg 1 ' 4 ) 


M 5 


Highest mass 


Mass gap 


Overclosure 


1.5- 1(T 6 


17 


40 keV 


10" 15 eV 


BBN 1 


3 ■ 10 13 


1.9 • 10 4 


1 GeV 


10 ^eV 


BBN 2 


1.7- 10 9 


1.0 • 10 5 


100 MeV 


1 meV 


BBN 3 


9 ■ 10 13 


9- 10 4 


1 GeV 


1 meV 


CMB 1 


4 • 10 10 


4- 10 4 


300 MeV 


0.1 meV 


CMB 2 


2 ■ 10 6 


1.3 • 10 4 


30 MeV 


1 fieV 


N(y) 


7- 10 s 


1.5 • 10 3 


100 MeV 


1 neV 


X-ray 


0.1 


7- 10 3 


500 keV 


0.5 fieV 


Reionisation 


2 


3- 10 4 


1 MeV 


10 ^eV 



Misalignment production 





Abundance (times M 21 ^ 4 ) 


M 5 (over 6^ 21 ) 


Highest mass 


Mass gap (over 24,/21 ) 


Overclosure 


1.3 • 10 27 


3.3- 10 7 


40 keV 


40 keV 


BBN 1 


2 ■ 10 37 


7- 10 8 


1 GeV 


300 MeV 


BBN 2 


1.2 ■ 10 35 


8 ■ 10 9 [5 • 10 8 ] 


100 MeV 


500 GeV 


BBN 3 


4 ■ 10 37 


3 • 10 9 [1.0 • 10 9 ] 


1 GeV 


30 GeV 


CMB 1 


7 ■ 10 35 


3 ■ 10 9 [7 • 10 8 ] 


300 MeV 


10 GeV 


CMB 2 


4 ■ 10 33 


1.9 • 10 9 [3 • 10 8 ] 


30 MeV 


10 GeV 


N{u) 


7 ■ 10 34 


1.4 ■ 10 8 


100 MeV 


3 MeV 


X-ray 


5 ■ 10 29 


5 ■ 10 9 [8 ■ 10 7 ] 


30 GeV 


100 keV 


Reionisation 


2 ■ 10 30 


1.1 • 10 10 [1.0 ■ 10 8 ] 


1 TeV 


500 keV 



Supernova 



Mb > 1.5 ■ 10 7 



TABLE I: Summary of all the constraints analysed in section IIVI The abundances and limits on M5 are expressed in TeV 
units. In square brackets we show the actual lower limit on M5 whenever the minimum mass gap demanded by the constraint 
is comparable to the highest mass considered in the same band. Indeed, for large initial O, the mass gap is pushed beyond the 
constraining band, whose upper limit corresponds to the 5D Planck mass shown in square brackets. The angle for which this 
happens can be found by comparing the last two columns. 



G. Astrophysics - Supernova bound 



The usual supernova bound arises |4| from the estimate of the energy loss by the star in the explosion, due to the 
escape of weakly interacting axions. The only difference is now that we have TR states (where T ss 30 MeV is the 
SN temperature) and the 4D limit (fixing /j = 10 7 TeV as before) becomes 



A 

/Tr 



> 10 5 TeV . 



(60) 



that is, M 5 > 1.5 • 10 6 TeV, valid as long as TR > 1, i.e. M 5 < 2 • 10 s TeV. 
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V. SUMMARY AND CONCLUSIONS 



In this work we have analysed the cosmology of 5D axions models, in particular in connection with brane world 
cosmologies. The main motivation for studying this problem comes from the fact that every time a new bulk field is 
introduced, the 4D cosmology is automatically affected by the appearance of new KK states; at the same time the 
high energy behaviour of the scale factor is known to be different than ordinary FRW cosmology, which has a profund 
impact on the abundances of particles generated during that epoch. 

In addressing this problem we have computed the abundances of KK axions produced by the three main (and most 
model independent) mechanisms: thermal production, the revitalised gravitational production, and production from 
misalignment. It has been shown why we expect that thermal abundances of KK modes is not likely to play any 
interesting role, whereas other production mechanisms can put stringent constraints on the parameters of the 5D 
model, of which we have chosen to single out the 5D Planck mass M5. 

Depending on some details of the 5D axion model and its precise (but unknown) cosmologically history around 
the inflationary epoch (if any), it is possible to constrain M5 to be not below some 10 4 TeV to 10 8 TeV, for a fixed 

= 10 7 TeV. This fiducial value for f± has been chosen since the usual 4D limits are most likely to hold in the 5D 
setup, especially in 5D, unless some strong fine tuning on the 5D parameters is required. 

Although there is no significant improvement as far as the 4D axion decay constant is concerned, we were still able 
to provide interesting lower limits on M5 (or /s, which is equivalent) that, for instance, preclude the existence of a 
TeV range (or even 10 3 TeV) extra dimension if the axion is allowed to propagate in the bulk. This conclusion relies 
heavily on the existence of a period of brane expansion in the early Universe, during which the abundances of KK 
particles can be greatly enhanced. The only way to completely evade the conditions we have reported is thus to build 
a model in which the Universe never had equilibrium energy densities high enough to allow for the aforementioned 
non standard expansion era. Of course even in this case some production mechanism is still available (gravitational 
particle production) but it is impossible to draw any conclusion, lacking a precise formulation of such cosmological 
epoch. Notice finally that in this paper we have not considered any solitonic production mechanism, which makes our 
limits conservative. 

Astrophysics is also known to be sensitive to the axion properties; if the axion field is comprised of a tower of 
different masses, and at least some of them contribute to supernova (or stellar) cooling, then it is possible to constrain 
M5 in yet another independent way. Using fi = 10 7 TeV we have found that supernova? limit M5 from below at 
about 10 7 TeV, which is, in some cases, the most relevant one. 

Before concluding, a note on the role of f± in this analysis. We have fixed this parameter before discussing the 
cosmological history of the KK axions, but it is easy to see what would have happened, had we picked up a different 
value for it. Consider first a larger f^. This means that axion interactions are weaker, and that the same KK 
state would decay later, rendering its contribution potentially more worrisome. This is indeed the case, and one can 
see that tighter limits would be obtained in this case (although if we confine ourselves to the 4D allowed range Q] 
the differences are mild). Had we chosen a much smaller /j instead the constraints coming from gravitational and 
misalignment productions would be relaxed. However, in that case thermal equilibrium would come back in, as the 
window at which the number of equilibrium degrees of freedom is dominated by Ukk broadens to much smaller M5, 
which means that (despite the faster decay rate) the enormous number of KK states reaching equilibration would 
survive well after Is, and the 5D Planck mass would need to be again pushed up enough for the first state to decay 
before BBN. 
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